Directional spatial structure of dissociated elongated molecular condensates 
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Ultra-cold clouds of dimeric molecules can dissociate into quantum mechanically correlated con- 
stituent atoms that are either bosons or fermions. We theoretically model the dissociation of cigar 
shaped molecular condensates, for which this difference manifests as complementary geometric struc- 
tures of the dissociated atoms. For atomic bosons beams form along the long axis of the molecular 
condensate. For atomic fermions beams form along the short axis. This directional beaming sim- 
plifies the measurement of correlations between the atoms through relative number squeezing. 

PACS numbers: 03.75.-b, 03.65.-w, 05.30.-d, 33.80.Gj 



The difference between bosons and fermions is funda- 
mental in physics. Ultra-cold degenerate quantum gases 
enable new kinds of explorations and applications of this 
difference. For example in Hanbury Brown and Twiss 
type measurements bosonic 4 He atoms are observed to 
be bunched while fermionic 3 He atoms are anti-bunched 
Besides such local spatial correlations, non-local cor- 
relations between pairs of atoms with opposite momenta 
have also been observed for both fermionic atoms 0] and 
bosonic atoms (3| . Such atom pairs, generated by disso- 
ciation of molecular-dimers [2|, or by four- wave mixing 
3 , are predicted to be entangled 0, 0, [f| 0, S, El EH EH • 
However the existence of non-classical Einstein- Podolsky- 
Rosen type correlations [4, 8] between atom pairs in these 
systems remains to be experimentally confirmed. 

We report an analysis of the effect of particle statistics 
on the dissociation of highly anisotropic Bose-Einstein 
condensates (BECs) of molecular dimers. When the dis- 
sociated atoms are bosons we confirm the expected beam- 
ing of atoms along the long axis of the cigar shaped 
molecular condensate by stimulated emission [7j. How- 
ever, when the atoms are fermions we find an unexpected 
beaming in the orthogonal direction, i.e., along the short 
condensate axis. The effect is due to Pauli blocking fol- 
lowed by atom-atom recombination which reduces the 
atomic density along the long axis. We also find enhanced 
relative number squeezing between the atoms in the half- 
spaces occupied by the oppositely directed beams. 

Directionality effects due to bosonic amplification have 
been observed [Til, EH in superradiant light scattering 
from elongated condensates @, EH EH- The fermionic 
counterpart of the effect due to Pauli blocking has not 
been discussed before to the best of our knowledge. 

The effective quantum field theory Hamiltonian de- 
scribing the system is given, in a rotating frame, by (IB ] 



H = H +h / d 
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Here, ^j(x,t) (i = 1,2) are the field operators for the 
atoms, which are two different spin states of the same iso- 
tope of mass m, and may be either bosonic or fermionic. 
The field operators satisfy the respective commuta- 
tion or anti-commutaion relations, [^(x, t), ^(x', t)] = 

<%5 2 (x-x') and {$ 4 (x,<),$](x',t)} = %<5 2 (x-x'). In 
the following we shall only consider two spatial dimen- 
sions, with x = (x,y). However we expect our results 
to be at least qualitatively valid for three-dimensional 
systems. 

The first term in Eq. ([T]) describes the kinetic energy 
of the atoms H = J dx£ i=1 2 ?i 2 | V$i| 2 /2m. The de- 
tuning A is defined so that spontaneous dissociation of 
molecules corresponds to A < 0, with 2fi|A| being the 
total dissociation energy that is converted into the ki- 
netic energy of atom pairs 17, 18j]. For molecules at rest, 
the dissociation primarily populates the resonant atomic 
modes in the two spin states having equal but opposite 
momenta, hh\ — — fik.2, with the absolute wavenumber 
equal to ko = |ki| = |k 2 | = y/2m\A\/h. 

In our analysis we use the undepleted molecular field 
approximation. It is valid for short dissociation times 
during which the converted fraction of molecules does 



not exceed about 10% [171 LUi ■ I n the regime of valid- 



(1) 



ity of the undepleted molecular approximation, the dis- 
sociation typically produces low density atomic clouds 
for which the s-wave scattering interactions are negli- 
gible figj | : hence their absence from our Hamiltonian. 
Additionally, the atom-molecule interactions can be ne- 
glected if the respective interaction energy per atom is 
much smaller than the total dissociation energy 2fi,|A|. 
We have confirmed the validity of our approximations 
(which improve with increasing |A|) in the case of dis- 
sociation into bosonic atoms by comparing the present 
results with exact, first- principles simulations using the 
positive- P representation [l9j]. For the case of fermionic 
atoms, first-principles simulations with multimode inho- 
mogeneous molecular condensates remain under develop- 
ment [20j |. and such a comparison is not presently possi- 



2 



ble. 

We assume that the bosonic molecules are in a coherent 
state initially, with the density profile Pm( x ) given by the 
ground state solution of the Gross-Pitaevskii equation 
with an anisotropic harmonic trap. We then have an 
effective, spatially dependent coupling 



fif(x) = xVpm(x) 



(2) 



The coupling coefficient x 13 is responsible for coherent 
conversion of molecules into atom pairs, e.g. via optical 
Raman transitions, an rf pulse, or a Feshbach resonance 
sweep 2l|, |22|, |23|, |2J, |25j. We assume that once the 
dissociation is switched on at time t = 0, the trapping 
potential is switched off, so that the evolution is tak- 
ing place in free space. Thus the role of the trapping 
potential is reduced to defining the initial shape of the 
molecular BEC. 

The Heisenberg equations for the atomic fields in the 
Hamiltonian JTJ) are then: 



a*i(x, t) 



dt 



A V 2 - a' 

2m 



$ 1 (x,t)± 5 (x)$^(x,i), 



(3) 



-^V 2 -A 

2 m 



$t( x , t)+ ff (x)*i(x, i). 



The +/— in the first equation correspond to 
bosonic/fcrmionic atoms. 

Expanding in plane wave modes, \¥j (x, t) = 
J G? 2 k aj (k, t) exp(— ik • x)/27r, where the operator ampli- 
tudes satisfy commutation or anticommutation relations, 
[ai(k,t),at(V,t)] = %<5 2 (k-k') or {a i (k,t),at(k',t)} = 
5ijS 2 (k— k'), according to the underlying statistics. The 
partial differential equations ([3]) reduce to a set of linear 
ordinary differential equations, 



' / " l(k ' /i -iA k S 1 (Kt)± J ^9(ci + k)aUq,t), (4) 
iA k al(k,t)+ / — 5(q-k)ai(-q,t), 



dt 

da^k, i) 



where <?(k) = J d 2 xg(x) exp(ik • x)/2tt is the Fourier 
transform of g(x) and = Hk 2 / (2m) + A, where k = 
|k|. For vacuum initial conditions the non-zero second- 
order moments are the normal and anomalous densi- 
ties; rij(k, k',i) = (aj(k, t)ai(k', t)) and mi 2 (k, k',t) = 
(ai(k, t)a2(k', f)). Higher-order moments can be ob- 
tained from these second-order moments using Wick's 
theorem, as the Hamiltonian is quadratic in the field op- 
erators in the undepleted molecular approximation. 

In a finite quantization volume the wave-vector k is 
discrete and the plane wave mode annihilation and cre- 
ation operators may be organized into a vector a. The 
Heisenberg equations ([4]) may then be written in vector- 
matrix form as da/dt = Ma, where M is a matrix of 
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Figure 1: (color online) Density in position space for the 
initial molecular condensate. It corresponds to the ground 
state of a harmonic trap with TF radii: Rtf,x — 60 ^im and 
Rtf, v — 6 [im. 
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Figure 2: (color online) Atomic density in momentum space 
(in units of m 2 ) at different times t after the start of dissocia- 
tion. The left and right columns are for bosonic and fermionic 
atoms, respectively. We have chosen the times ti/to = 1.1 for 
(a) and (b), ti/to = 2.2 for (c) and (d), and tz/to = 3.3 for 
(e) and (f). Here to = 1 / (Xy/ Pm.o) = 2.5 ms is the time scale. 



complex numbers. Linearity ensures that the solutions 
of these operator equations can be found by numerically 
computing the matrix exponential exp(Mt). We used an 
81x81 numerical lattice in momentum space, which gave 
the same results as for a 61 x 61 lattice. 

Here we explore the role of a strongly anisotropic sys- 
tem by assuming that the initial molecular condensate 
was trapped in a harmonic potential with a frequency 
along the y axis u) y that is ten times that along the x axis; 
Lo y = lQuj x - We assume a molecule-molecule scattering 
length such that the Thomas-Fermi (TF) approximation 
is valid. The corresponding molecular density profile is 
Pm = Pm,q (l — x2 /Rtf,x 



y 2 /R 2 :F y), where the Rtf, 
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arc the TF-radii in each direction. The parameter values 
we used are given in Ref. (26j . 

Fig. [5] shows the density of our initial two-dimensional 
molecular condensate in position space. Fig. [2] shows 
the total density of the dissociated atoms in momentum 
space, at three times as dissociation progresses. The left 
column shows bosonic atoms, and the right column shows 
fermionic atoms. Note that after sufficient time-of-flight 
expansion, these momentum s pac e distributions are re- 
produced in position space 0,113]. At the earliest time 
there is little difference between the bosons and fermions, 
and the width of the distribution is determined by the 
energy-time uncertainty relation. Later, the spatial dis- 
tributions of the bosons and fermions develop quite dif- 
ferently. For bosons, our results can be compared with 
first-principles simulations using the positive- P represen- 
tation (lij , in which the molecular field and its depletion 
is treatedquantum mechanically; the comparison is given 
in Ref. [27j and shows good agreement. 

The highest densities of bosons develop along the long 
axis of the molecular condensate. This is due to Bose 
stimulation, which in the undepleted molecular field ap- 
proximation leads to approximately exponential growth 
of the atom number. In contrast, the highest densities of 
fermions occur in the orthogonal direction; that is along 
the short axis of the initial molecular condensate. Phys- 
ically, the effect is due to Pauli blocking followed by the 
reduction of the atomic density due to atom-atom re- 
combination along the long axis. The dynamics along 
the short axis (for sufficiently small Rtf.v), on the other 
hand, does not reach the regime dominated by satura- 
tion of mode populations and atom-atom recombination 
as the atoms propagating along y leave the molecular 
condensate earlier than those propagating along x. 

Mathematically, the difference is due to the sign differ- 
ence on the right-hand-side of the Heisenberg equations 
([3|), corresponding to bosons and fermions. As noted 
in Ref. [28| the sign difference determines sinh(oifcT) 
population growth for bosons, and sin(afcT) popula- 
tion oscillations for fermions; where t = t/to is a di- 
mensionless time, with to — 1/{X\/Pm,o), and a>k — 
\Jl ± t\{hk 2 1 (2m) + A) 2 ; here, the + corresponding to 
fermions and the — to bosons. 

The fermionic oscillations - both in momentum space 
and in time - can be seen in Figs. [2] (d) and (f) and 
in Fig. [31 As the solutions of Ref. [28| are for spatially 
homogenous systems they do not quantitatively describe 
these minima for our inhomogeneous molecular conden- 
sate. Nevertheless they predict aspects of the qualitative 
behavior, such as the movement of the density minima 
and maxima to lower values of k x with increasing time; 
Figs. [2] (d) and (f). The single density maximum along 
the k y axis is a result of the small TF radius along that 
axis, Rtf,v — 6 /im. The atomic velocity at |/c/fco| = 1 
is 0.95 mm s _1 , so during t\ = l.ltp = 1.1(2.5 ms) the 




Figure 3: (color online) Atomic density in momentum space 
at the position of resonance (|k| = ko) along the two dif- 
ferent Cartesian directions, in units of 10 -11 m 2 . The two 
upper (red) curves are for bosons which grow exponentially 
due to Bose stimulation. The two lower (blue) curves are for 
fermions. Solid lines are along the A^-direction (k x = ko, k y = 
0), dashed lines are along the fc y -direction (k x — 0,k y — ko). 
The inset shows the fractional atom number (relative to the 
total initial number of molecules), in units of tenths of a per- 
cent, for fermions (solid blue) and bosons (dashed red). 

atoms travel 2.6 fim, and during = 3.3to they travel 7.8 
/xm, which is greater than the TF radius Rtf.v- Hence 
a particular atom traveling in the y direction typically 
interacts with the molecular condensate for less than the 
time £3. This accounts for the similarity of Figs. [2] (d) 
and (f) along the y axis. 

We have used our solutions of the Heisenberg equa- 
tions (TJ|) to calculate the correlations between atom num- 
ber fluctuations in the two different spin states 1 and 2. 
We considered momentum areas centered on the opposite 
resonant momenta, fck.Q and — ftko, for the smallest and 
largest possible areas in which the atoms are captured 
and their number is measured. The smallest momentum 
areas are the numerical lattice areas, and the largest are 
opposite halves of the momentum space. These are ori- 
entated differently for fermions and bosons, guided by 
Figs. H](e) and (f), to capture one beam in each half: for 
bosons the halves are split by the k y axis and for fermions 
by the k x axis. 

For the case of the smallest momentum area (the 
area of the numerical lattice cell, Ak x Ak y ), we de- 
fine the atom number operators via iV/,±ko (t) = 
?ij(±ko, t)Ak x Ak y , where ko = k$e x for bosons and 
ko = kge y for fermions, with e x and e y being the Carte- 
sian unit vectors. We quantify the correlations by the 
normalized relative atom number variance, 

Vko.-kc (*) = ([A(n 1)ko - n 2 ,_ ko )] 2 )/A S N, (5) 

where AC = C — (C) is the fluctuation in C and Asn = 
((Ani i k ) 2 ) + ((An2,-k ) 2 ) defines the uncorrelated shot- 
noise level. Variance Vk ,-k (t) < 1 implies squeezing 
of the relative number fluctuations below the shot-noise 
level. The results (see Ref. 27] for details) are shown 
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produce qualitatively different geometrical distributions 
for bosonic and fermionic atoms. These are in the form 
of "twin" beams in opposite half-spaces. The squeezing 
of the relative atom number fluctuations between these 
beams can be quite strong. 

The authors acknowledge support from the Australian 
Research Council through the ARC Centre of Excellence 
scheme. MO acknowledges IPRS/UQILAS. 
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Figure 4: (color online) Relative number variance according 
to Eqs. ([5jl and (f6]l. raw (unbinned) and binned, for fermions 
(solid blue) and bosons (dashed red) . The dotted (black) line 
in the inset is the t <C to asymptote, Eq. @ . 



in Fig. [5] (upper curves), where we see relatively small 
degree of squeezing for this "raw" (unbinned) case. 

The squeezing in the relative number can be enhanced 
for larger counting areas, which we refer to as binning 
[Io| . For the largest pair of bins coinciding with the 
left (L) and right (R) halves of the momentum space 
for bosons, and with the bottom (B) and top (T) halves 
for fermions, we can introduce particle number operators 



N 



L(R) 



and iV?^ {i — 1,2). The normalized variance of 
the relative number fluctuations between the spin states 
1 and 2 is defined similarly to Eq. ([5]), 



Vt 



LR(BT) 



(t) = ([A(J\£ 



CrR(T) 



At (B) )] 2 )/A 



SN- 



(6) 



The results for Vlr(bt)W ( see Ref. I27J for details) are 
shown in Fig. 21 comparing the upper and lower curves, 
we see that binning into momentum half-spaces reduces 
the relative number variance, or increases the degree of 
squeezing, by a factor of about ten for fermions and a 
hundred for bosons. This difference is due to the different 
momentum uncertainties along the long and short spa- 
tial axes. The momentum uncertainty is about ten times 
smaller along the long axis, which is the bosonic beaming 
direction. Hence there are about ten times less correlated 
pairs ending up in the same momentum half-space. For 
the fermions, on the other hand, the relatively high mo- 
mentum uncertainty in the beaming direction produces 
more pairs in the same half-space, reducing the number 
difference squeezing. 

In the case of unbinned atom numbers the fermionic 
and bosonic cases have a common short time asymptote 
that can be determined analytically using a perturbative 
approach [29j. Applying it to our 2D system gives 

Vk ,- ko = 1 - 2R T F,xRTF, y Ak x Ak y / (9tt) ~ 0.978, (7) 

which agrees with our numerical results (see Fig. 2]). 

In summary, we have shown that dissociation of elon- 
gated condensates of molecular dimers into atoms can 
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1. Comparison with the positive-P simulations. 
Here, we compare the results of the present treatment 
using the undepleted molecular approximation and exact 
first-principles simulations using the positive-P represen- 
tation 2, 1|| , in which the molecular depletion is taken 
into account within the full quantum field theory Hamil- 
tonian. In Fig. [5] we show the corresponding atomic den- 
sities for the case of bosonic atoms. Figure [5] (a) is the 
same one as Fig. 2 (e) of the main text, while Fig. [5] (b) 
is the corresponding positive-P result. As we see the 
agreement is very good. First principles simulations for 
fermionic atoms using the Gaussian stochastic methods 
of Ref. [2(| are still under development, but we expect a 
similarly good agreement as the_ molecular depletion is a 
weaker effect in this case 17 
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Figure 5: Atomic density in momentum space (in units of 
m 2 ) at ts/to =3.3 within the undepleted molecular field ap- 
proximation (a) and using the positive-P simulations with 
molecular depletion (b). 



2. Shot noise and relative number variances. — For 
bosons the uncorrelated shot-noise is determined by Pois- 
sonian statistics so that in the case of the smallest count- 
ing area we have 

A SN = ((Anx, ko ) 2 ) + <(An 2 ^ k0 ) 2 ) 

= <«i,k ) + (n2,-k D >- (8) 

For fermions, on the other hand, the shot-noise is always 
sub-Poissonian and is given by [28| 

A SN = ((An^J 2 ) + <(An 2 ,_k ) 2 > 

= (ftl,ko>(l - (*H,ko)) + ("2,-0(1 - (»2,-ko»- 

(9) 

The relative number variance can be written in the 



t 7 ,.s ("i,k n "i,kj ~ (n 1M n 2 ,-k ) rin , 
Kk ,-k (tj = t , (1U) 

where we have taken into account that (rii. ko ) = 
("2,-ko) = «k and (ni,k ni, ko ) = (n 2 ,-k n 2 ,- ko ). Note 
that (ni iko ni jko ) = (fta jko ) = n ko for fermions. Ap- 
plying Wick's theorem to calculate the higher-order mo- 
ments, and introducing the resonant anomalous density 
™k = (ai,k a2,-k )j where a ijk (t) = %(k, t)yfAk x Ak y 
we obtain 



V ko ,- ko (t) = l- 
for bosons, and 

V ko ,_ ko (t) = l- 



l^kol 



"kn 



l TO k | 



n ko (l-n ko ) 



(11) 



(12) 



for fermions. The variance V koj _ ko (t) is then calculated 
by numerically solving Eqs. (4) of the main text. 

For the case of the largest pair of counting areas, which 
we refer to as binning, we consider the left (L) and the 
right (R) halves of the momentum space for bosons, and 
the bottom (B) and top (T) halves for fermions. Ac- 
cordingly, we introduce particle number operators jvH R ^ 
and Nf [T) (j = 1,2), defined as A> L = £ L h jM (with 
a similar convention for R, B, and T), where X^ L and 
stand for double sums = Sfc x <oSfc an d 

Y, kx T,k y <o and 
The normalized variance of the rel- 



E R = Efc x >oEfc„> whereas £ B 



ative atom number fluctuations in the spin states 1 and 
2 is defined in Eq. (6) of the main text, in which the 
uncorrelated shot-noise level 



A SN - ^((Anr.k) 2 ) + ^((An 2 , k ) 2 ) 



(13) 



for bosons is given by 

A S n = £ R <^) + = <^i R ) + <^ L >, (14) 

while for fermions it is given by 

A S n = ^ T (ni, k )(l - (n x , k )) + ^ B (n 2 , k )(l - (ri 2 ,k», 

(15) 

which we note is not same as [(A\ T ) (1 - (A\ T ) ) + (iVf 
(N®) )] • Applying Wick's theorem to factorize the higher- 
order moments in Eq. (6) of the main text, we obtain 



Z! k , 1 ( a lk a i.k') i 2 - Y?L i ( a i.k a 2,k 



(16) 
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for bosons, and 

VkT(t) 



" Yll I ( a lkOl,k'> I 2 - Yll Y k , I (0l,k02,k'> i 



for fermions. Here, we have used the fact that (fti,k) — 
(»i,-k) = (n 2 ,k), (N?) = (iV 2 L ), and (iVf) = (7V 2 B )', due 
(iV] r ) (17) to symmetry considerations. The calculation of the vari- 

ances Vlr(£) and Vbt(£) can now be done using the nu- 
merical solution of Eqs. (4) of the main text. 



